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ABSTRACT 

The dense cores of X-ray emitting gaseous halos of large elliptical galaxies with temperatures 
kT < 0.8 keV show two prominent Fe xvn emission features, which provide a sensitive diag- 
nostic tool to measure the effects of resonant scattering. We present here high-resolution spec- 
tra of five bright nearby elliptical galaxies, obtained with the Reflection Grating Spectrometers 
(RGS) on the XMM-Newton satellite. The spectra for the cores of four of the galaxies show 
the Fexvii line at 15.01 A being suppressed by resonant scattering. The data for NGC 4636 
in particular allow the effects of resonant scattering to be studied in detail and to prove that 
the 15.01 A line is suppressed only in the dense core and not in the surrounding regions. Us- 
ing deprojected density and temperature profiles for this galaxy obtained with the Chandra 
satellite, we model the radial intensity profiles of the strongest resonance lines, accounting for 
the effects of resonant scattering, for different values of the characteristic turbulent velocity. 
Comparing the model to the data, we find that the isotropic turbulent velocities on spatial 
scales smaller than ^1 kpc are less than 100 km s"' and the turbulent pressure support in the 
galaxy core is smaller than 5% of the thermal pressure at the 90% confidence level, and less 
than 20% at 95% confidence. Neglecting the effects of resonant scattering in spectral fitting 
of the inner 2 kpc core of NGC 4636 will lead to underestimates of the chemical abundances 
ofFe and O by -10-20%. 

Key words: X-rays: galaxies - X-rays: individual: NGC 4636, NGC 4649, NGC 5813, 
NGC 1404, NGC 4472 - galaxies: abundances - cooling flows - turbulence 



1 INTRODUCTION 

The hot intra-cluster medium and the hot halos around giant ellipti- 
cal galaxies are usually assumed to be optically thin. Although this 
assumption is valid for most of the emitted X-ray photons, at the 
energies of the strongest resonant transitions, the hot plasma can 
be optically thick (Gilfanov et al. 1987). The transition probabili- 
ties of strong resonance lines are large and, if the column density 
of the ion along a line of sight is sufficiently high, photons with 
the energies of these resonance lines will get absorbed and, within 
a short time interval, reemited in a different direction. Because of 
the short time between absorption and emission, this process can 
be regarded as scattering. Since resonant scattering in clusters and 
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elliptical galaxies will cause the radial intensity profile of an emis- 
sion line to become weaker in the centre and stronger outside, it 
can also lead to an underestimate of metal abundances in the dense 
cores and a corresponding overestimate in the surrounding region. 
This e ffect is, however, not expected to be large. , Sanders & Fabiai] 
( l2006h found that metallicities in cluster cores could be underesti- 
mated by at most 10% due to resonant scattering. 

iGilfanov et al. ||l983) pointed out that since the optical depth t 
in the core of a resonance line depends on the characteristic velocity 
of small-scale motion, measurements of r give important informa- 
tion about the turbulent velocities in the hot plasma. Measuring the 
level of resonant scattering in clusters and giant elliptical galaxies 
is thus a good way to constrain the energy in turbulence and the tur- 
bulent pressure support. The first constraints on turbulent velocities 
using resonant scattering were obtained for the Perseus cluster by 
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Figure 1. Fraction of Fe in the form of Fe xvn and the theoretical line ratio 
(^il7.05 + ^il7.io)/Au5.oi for an optically thin plasma as a function of the 
plasma temperature. The vertical lines mark the RGS measured tempera- 
tures of the five galaxies in our sample. 

IChurazovet"aLl ( l2004h . Using XMM-Newton EPIC data, they com- 
pared the relative fluxes of the ls-2p and ls-3p He-like Fe lines in 
the core and in an annulus surrounding the core of the Perseus clus- 
ter. The expected optical depth of the 6.7 keV ls-2p Fexxv reso- 
nance line is much larger than that of the ls-3p line, therefore the 
ratio provides inform ation about the level of resonant scattering. 
IChurazovetaLl ( l2004h found no evidence for resonant scattering in 
Perseus, indicating that differential gas motions on scales smaller 
than ~100 kpc in the core of the cluster must have a range of ve- 
locities of at least half of the sound speed. Independently, using the 
same data, iGastaldello & Molendil (j2004) reached similar conclu- 
sions. 

The most sensitive spectral lines to determine the level of res- 
onant scattering are the Fexvn lines at 15.01 A (2p-3d) and the 
unresolved blend of the same ion at 17.05 and 17.10 A (2p-3s). 
The optical depth is directly proportional to the column density of 
the ion and the oscillator strength of the given transition. While the 
oscillator strength of the 15.01 A line is / = 2.73, and thus the line 
is expected to have a relatively large optical depth, the oscillator 
strength of the 17.05 A line is / = 0.12 and the optical depth of 
this line blend is negligible (the oscillator strength of the 17.10 A 
line is of the order of lO"**!]. Because of this dramatic difference 
in expected optical depths, and the fact that both lines originate 
from the same ion of the same element, the Fexvn lines provide, 
in principle, an excellent diagnostic tool to measure the magni- 
tude of resonant scattering. At the relatively low temperatures of 
hot halos around giant elliptical galaxies, ~0. 6-0.8 keV, both lines 
are very strong and, in this temperature range, their expected in- 
tensity rati os have only a very w eak dependence on temperature 
(see Fig.m'P oron & Behaill2002h . Therefore, any observed spatial 
dependence of the intensity ratios is unlikely to be due to gradi- 
ents in temperature or Fe abundance in the hot plasma, making 
these lines especially clean diagnostic tools. The Fexvii line ra- 
tios in extended objects can be measure d with the XMM-Newto n 
Reflection Grating Spectrometers (RGS, Iden Herder et ^ l200lh . 

' http://cxc.harvard.edu/atomdb/WebGUIDE/index.html 




Figure 2. Chandra image of NGC 4636 with the RGS extraction regions 
over-plotted. The central extraction region is 0.5' wide (indicated by full 
lines), the two extraction regions next to the core (between the full and 
dashed lines) are 2.25' wide. The image was extracted in the 0.5-2.0 keV 
band, cleaned and adaptively smoothed. 



IXu et al.l ( l2002h reported the first and so far only radial profile of 
the (/ii7.05 + /.ii7.io)//,u5.oi ratio for the X-ray luminous elliptical 
galaxy NGC 4636, which indicates the presence of resonant scat- 
tering. Attempts to measure resonant scattering in NGC 5044 and 
M 87, which are hotter and mo re massive systems with weaker 
Fe xvn lines, were not successful jTamura et alj2003l : IWemer et al.l 
2006). 

Here we follow up on the work bv lXu et aT H200i), by report- 
ing the Fexvn (/^n.os + lAn.K))lhi50\ line ratios measured in the 
cores of five nearby bright elliptical galaxies observed with XMM- 
Newton RGS: NGC 4636, NGC 5813, NGC 1404, NGC 4649, and 
NGC 4472. In Sect. 2 we describe the sample and observations; in 
Sect. 3 we discuss the details of the data analysis; and in Sect. 4 
we present the results of the observations. We focus our attention 
on NGC 4636, for which data are of very high quality. Because 
NGC 4636 is X-ray bright, has a relatively long exposure and a 
favorable temperature, for this system we can measure and com- 
pare the observed line ratios in the core and surrounding regions. In 
Sect. 5 we present a model for the radial profile of the 15.01 A line 
in NGC 4636 for different values of turbulent velocities, derived 
using deprojected density and temperature profiles from Chandra 
data. Comparing this model to the data, we place constraints on 
the characteristic velocity of isotropic turbulence in the core of the 
galaxy. Implications of the results are discussed in Sect. 6. 

Throughout the paper, abundances ar e given with respect to 
the "proto-solar values" by iLoddersI |_2003). All errors are quoted 
at the 68% confidence level for one interesting parameter (AC = 1; 
for C statistics). 



2 SAMPLE AND XMM-NEWTON OBSERVATIONS 

We have analyzed a sample of five nearby, giant elliptical galaxies 
observed with XMM-Newton RGS. The sample is listed in Table[T] 
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Table 1. List of the galaxies in the sample, with their distance (luminosity distances, o btained from the N ASA/IPAC Extragalactic Database, 
|http://nedwww.ipac.caltech.edu/} , corresponding linear scale per arcminute. Galactic Wh value jKalberla et alj EOO^ . average temperature measured by 
ROSAT. The last three columns give the XMM-Newton revolution of the observation, and the total and filtered RGS exposure times, respectively. 



galaxy 


distance 


scale 




kT 


rev. 


Exp. time 


Clean time 




Mpc 


kpc/arcmin 


lO^o cm-2 


keV 




ks 


ks 


NGC 4636 


17.5 


5.1 


1.90 


0.55* 


197 


64406 


58205 


NGC5813 


29.9 


8.7 


4.37 


0.52+ 


1029 


36708 


29899 


NGC 1404 


25.5 


7.4 


1.51 


0.60* 


1033 


55032 


17231 


NGC 4649 


19.8 


5.8 


2.1 


0.78* 


196 


54210 


46096 


NGC 4472 


18.4 


5.4 


1.53 


0.88* 


744 


111033 


81850 



• O'SuUivan et al. (20031 

t Reiprich & Bohringer (2002) 





Figure 3. Chandra images of the giant elliptical galaxies in our sample with over-plotted RGS extraction regions (see also Fig.|2). The extraction regions are 
0.5' wide. While the core of NGC 5813 is strongly disturbed, the cores of the other three galaxies are relatively relaxed. The images were extracted in the 
0.5-2.0 keV band, cleaned and adaptively smoothed. 



The galaxies were selected based on their relatively low core tem- 
peratures, at which a significant fraction of Fe is expected to be 
in the form of Fexvii (see Fig. [T}, their large X-ray fluxes, and 
strongly peaked surface brightness distributions, which are nec- 
essary to obtain RGS spectra with sufficient statistics. The sam- 
ple contains both relatively relaxed galaxies and systems showing 
very clear signs of interaction between the hot halo gas and cen- 



tral Active Galactic Nucleus (AGN). In Figs. |2] and [3] we show 
0.5-2.0 keV Chandra images of the galaxies in our sample with 
the RGS extraction regions overplotted. The images show that the 
cores of NGC 4649, NGC 4472, and NGC 1404 are relatively re- 
laxed (although on larger scales NGC 1404 shows a prominent 
cold front to the Northwest, and NGC 4649 and NGC 4472 harbor 
central radio sources, which exhibit interaction with the hot gas, 
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IShurkin et al.ll2008l : lBiller et alj2004l) . NGC 4636 and NGC 5813, 
on the other hand, look much more disturbed, with obvious signs 
of recent episodes of interaction between the hot plasma and the 
AGN. NGC 4636 has a very dense core and X-ray brig ht arm-like 
structures which might have been produced by shocks ( J ones et al] 
[2002). The core of NGC 5813 is highly disturbed with most of the 
hot gas being pushed into dense sheets by the rising bubbles of rel- 
ativistic plasma. 

Currently the best data set for the study of resonant scattering 
is that for NGC 4636. The X-ray halo of this galaxy has a conve- 
niently low temperature, at which ~26% of Fe is in the form of 
Fe xvii, a very dense core and a deep XMM-Newton observation, 
which also allows us to extract spectra from outside the core region 
with reasonably good statistics in the Fexvii lines. Using a deep 
Chandra observation of this galaxy we derive deprojected density 
and temperature profiles which we use to model radial intensity 
profiles for the strongest resonance lines, taking into account the 
effects of resonant scattering for different assumed values of turbu- 
lent velocities. 



3 DATA ANALYSIS 

3.1 XMM-Newton RGS data analysis 

We processed the RGS data using version 8.0.0 of the XMM- 
Newton Science Analysis Syst em (SAS), and extrac ted spectra us- 
ing the method described by iTamura et al. I ( l200lh . To minimize 
the contamination by soft-protons from Solar flares, we extracted 
a light curve for each dataset using events on CCD 9 of the RGS, 
outside the central region, with a distance larger than 30" from the 
dispersion axis, and excluded time intervals with elevated count 
rates. The original exposure time and the total time after cleaning 
are listed in Table [T] The flares in the observation of NGC 1404 
were weak and we obtain the same best fit spectral parameters for 
both the filtered and unfiltered data set. Therefore, for this galaxy 
we used the full observation. Because emission from the galaxies 
fills the entire field of view of the RGS, the background cannot 
be estimated from a region on the detector away from the source. 
Therefore, we modeled the background using the standard RGS 
background model for extended sources available in SAS (rgsbkg- 
model,|Gonzalez-Riestra 2004). 

The effects of resonant scattering are most significant in the 
centermost parts of the galaxies, where the density is highest. 
Therefore, we extract spectra from narrow 30" wide bands centred 
on the cores of the galaxies, probing regions within radii of 1.3- 

2.2 kpc (depending on distance) in the cross dispersion direction 
of the RGS. We extract and fit both first and second order spec- 
tra. Because the RGS operates without a slit, it collects all photons 
from within the 0.5' x ~12' field of view. Line photons originating 
at angle A6 (in arcminutes) along the dispersion direction, will be 
shifted in wavelength by 

= 0.138A6' A. (1) 

Therefore, every line will be broadened by the spatial extent of the 
source. To account for this broadening in our spectral model, for 
each RGS spectrum, we produce a predicted line spread function 
(LSF) by convolving the RGS response with the surface brightness 
profile of the galaxy derived from the EPIC/MOSl image in the 
0.8-1.4 keV band along the dispersion direction. Because the radial 
profile of a particular spectral line can be dilferent from the overall 
radial surface brightness profile (e.g. in the presence of metallicity 



gradients), the line profile is multiplied by a scale factor s, which is 
the ratio of the observed LSF width to the expected LSF for a flat 
abundance distribution. This scale factor is a free parameter in the 
spectral fit. 

For the spectral modeling of the RGS data we use the SPEX 
package ( Kaastra et al. 1996). Spectral fitting of the RGS data is 
done in the 10 A to 28 A band. We model the emission of the 
observed galaxies as optically thin plasmas in coUisional ioniza- 
tion equilibrium, absorbed by neutral Galactic gas with Solar abun- 
dances. For each object we fix the Galactic column density to the 
value deter mined by the Leiden/ A rgenitine/Bonn (LAB) Survey of 
Galactic H i jKalberla et alj2005l) . For the emission of the elliptical 
galaxies we assume a single-temperature plasma. Using more com- 
plicated multi-temperature models does not improve the fits and 
differential emission measure models always converge to a sim- 
ple single-temperature approximation. The spectral normalization, 
plasma temperature, and the abundances of N, O, Ne, and Fe are 
free parameters in the fit. For the spectral fitting we use C-statistics. 

3.2 Chandra analysis of NGC 4636 

We reprocessed the Chandra observations of NGC 4636 (OBSID 
3926 and 4415) applying the latest CTI and time-depen dent gain 
calibrations using standard methods (see lVikhlinin et al]|20 05. for 
more details). In brief, we performed the usual filtering by grade, 
excluded bad/hot pixels and columns, removed cosmic ray 'after- 
glows', and applied the VF mode filtering. We excluded data with 
anomalously high background. The remaining exposure time for 
NGC 4636 is 150.4 ks. The background files were processed in the 
same manner as the observations. An additional background com- 
ponent is produced during the 41 ms of the 3.2 s nominal exposure 
readout of the ACIS CCDs, while the chips are exposed to the sky. 
This small contribution of the source flux, 1.3%, is uniformly re- 
distributed along the re adout direction an d is sub tracted using the 
technique described bv iMarkevitch et ^ ^20001) . The blank-field 
background is also renormalized by reducing its integration time 

by 1.3% to account for this additional subtra ction. 

As the final step in data preparation (see lChurazov et alj2008l 
for more details), we computed, for each X-ray event, the ratio of 
the effective area (a function of photon energy and position) to that 
at a predefined energy and position: 

ri = A(E,xj,yd)/Ao{E), (2) 

where A(E, Xd,yd) includes mirror and detector efficiencies (includ- 
ing non-uniformity of the detector quantum efficiency and the time 
and spatially dependent contamination on the optical blocking fil- 
ter). Finally, for each event list, we make an exposure map that 
accounts for all position dependent, but energy independent, effi- 
ciency variations across the focal plane (e.g., overall chip geometry, 
dead pixels or rows, and variation of telescope pointing direction). 
We use these data in Sect.|5]to derive deprojected gas density and 
temperature profiles. 



4 OBSERVATIONS OF RESONANT SCATTERING 

The results from the spectral fitting of the XMM-Newton RGS data 
with an optically thin single-temperature plasma model are shown 
in Tableland Figures |4] and [5] The 13.8-15.5 A part of the spec- 
trum, where the strongest Fexvii and Fexvm resonance lines are 
present, was initially excluded from this fit. These lines are ex- 
pected to be suppressed by resonant scattering and this suppres- 
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Table 2. The best fit parameters for a single-temperature, optically thin plasma model fitted to the XMM-Newton RGS spectra extracted from 0.5' wide regions 
centred on the cores of the galaxies. For NGC 4636 we also show the results from fits to spectra extracted from two 2.25' wide regions surrounding the core 
(see Fig. |2). The 13.8-15.5 A part of the spectrum, where the strongest Fexvn and Fexvm resonance lines are present was initially excluded from the fits. 
Fluxes are given in the 0.3-2.0 keV band. The emission measure is defined as 7 = jnciiiidV. The scale factor s is the rati o of the observed LSF width to the 
expected LSF for a flat abundance distribution. Abundances are quoted with respect to the proto-solar values of lLodder3 i2003). The last three rows list the 
best fit line ratios in the full spectral band (after the Fe xvn ion was set to zero in the model and replaced by gaussian lines), the theoretical line ratios predicted 
for an optically thin plasma, and the derived level of suppression of the 15.01 A line, (///o)i5 qia- 



galaxy 


NGC 4636 core 


NGC 4636 outer reg. 


NGC 5813 


NGC 1404 


NGC 4649 


NGC 4472 


flux (10"'2 (.n,-2-, 


1.75 ±0.08 


2.57 ±0.17 


1.47 ±0.12 


1.08 ±0.11 


1.63 ±0.10 


1.44 ±0.08 




0.47 ± 0.02 


0.46 ± 0.03 


1.01 ±0.10 


0.50 ± 0.05 


0.31 ±0.03 


0.34 ± 0.02 


kT (keV) 


0.606 ± 0.006 


0.695 ± 0.004 


0.645 ± 0.008 


0.608 ± 0.009 


0.774 ± 0.007 


0.781 ±0.006 


s 


0.40 ± 0.04 


1.02 ±0.04 


0.87 ±0.11 


0.97 ± 0.22 


0.69 ± 0.21 


0.79 ±0.12 


N 


1.3 ±0.3 


1.5 ±0.4 


2.0 ±0.8 


2.3 ±0.8 


1.3 ±0.7 


1.3 ±0.5 





0.44 ± 0.05 


0.61 ± 0.06 


0.53 ± 0.09 


0.58 ±0.10 


0.61 ±0.15 


0.53 ± 0.07 


Ne 


0.31 ±0.08 


0.39 ±0.18 


0.33 ±0.19 


0.81 ±0.22 


1.31 ±0.35 


1.18 ±0.22 


Fe 


0.52 ± 0.03 


0.92 ± 0.06 


0.75 ± 0.09 


0.67 ± 0.08 


0.87 ±0.18 


0.83 ± 0.08 


[(.^,117.05 + ^,117.10)/^il5.0l]obscrved 


2.04 ± 0.21 


1.28 ±0.13 


1.99 ±0.34 


1.98 ±0.29 


1.25 ±0.28 


2.24 ± 0.34 


[(^.117.05 + /,117.10)/^il5.0l]predictcd 


1.31 


1.31 


1.30 


1.31 


1.27 


1.27 


hsm/Io 15.01 


0.64 ± 0.07 


1.02 ±0.10 


0.65 ±0.11 


0.66 ±0.10 


1.02 ±0.23 


0.57 ± 0.09 




14 16 
wavelength (A) 



14 16 
wavelength (A) 



Figure 4. RGS spectra for NGC 5813, NGC 1404, NGC 4649, and NGC 4472. The spectra were extracted from 0.5' wide extraction regions centred on the 
cores of the galaxies. The full line indicates the best-fit optically-thin single-temperature plasma model. The 13.8-15.5 A part of the spectrum, where the 
strongest Fe xvn and Fe xviii resonance lines are present, was excluded from these spectral fits. 
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wavelength (A) wavelength (A) 



Figure 5. XMM-Newton RGS spectra extracted from a 0.5' wide extraction region centered on tlie core of NGC 4636 (left panel) and from two 2.25' wide 
extraction regions surrounding the core (right panel). The ratios of the Fe xvn lines at 15.01 and ~ 17 A change. The full line indicates the best fit optically thin 
single-temperature plasma model. The 13.8-15.5 A part of the spectrum, where the strongest Fexvn and Fexviii resonance lines are present, was excluded 
from these spectral fits. 



sion, which is not accounted for in the plasma model, could slightly 
bias the best fit spectral parameters. After obtaining the best fit, we 
freeze the temperature and Fe abundance of the thermal component, 
set the abundance of the Fe xvii ion to zero in the model, add narrow 
Gaussians at the wavelength of the strongest Fe xvii lines (restframe 
wavelengths of 15.01 A, 17.077 A, 15.30 A, and 16.80 A) and re- 
fit the full 10-28 A spectral band. Freezing the temperature and 
the Fe abundance is necessary because by setting the Fe xvn ion to 
zero and replacing it with Gaussians in the model we lose important 
constraints on these parameters. The line normalisations, the con- 
tinuum level (normalisation of the thermal component), the abun- 
dances of the other elements and the broadening factor s are free 
parameters. In the error calculation, we also marginalise over the 
uncertainties of the temperature and Fe abundance as determined in 
the initial fit. We determine the best fitting normalisations for these 
Gaussians and the value of (/,u7.05 + /,ii7.io)//,ii5.oi (the 17.05 and 
17.10 A lines are blended and we fit them with a single Gaussian 
at 17.077 A). We divide the theoretical ratio (/.m os + /,ii7.io)/^ii5.oi 
by the measured value, which, assuming the lines at ~17 A are 
optically thin, gives us the level at which the 15.01 A line is sup- 
pressed, (///o)i5oiA- These values are listed in Table|2] In four out 
of five galaxies the line ratios are significantly higher (observed val- 
ues ~ 2.0) than the expected theoretical value for an optically thin 
plasma of 1.27-1.31, indicating that the 15.01 A line is optically 
thick. In NGC 4649, the line ratios are consistent with the 15.01 A 
line being optically thin. 

A number of factors argue that the observed high (/,i 17.05 + 
^,ii7.io)/^ii5.oi line ratio is due to resonant scattering. The (/jn.os + 
/,ii7.io)/^ii5.oi line ratio in an optically thin plasma can reach val- 
ues ~ 1.9 only for t emperatures as low as 0.2 keV (see Fig. [T] 
iDoron & Beharl2002l) . which is much lower than the observed tem- 
peratures of the galaxies. Moreover, if the plasma in the cores of 
the galaxies were multiphase with an additional 0.2 keV compo- 
nent, which would alter the Fe xvii line ratios, then we would also 
observe other strong lines (e.g. Ovii), which we do not. Another 
possibility to obtain such high line ratios in an optically thin plasma 
is non-equilibrium ionization (NEI). However, NEI effects are un- 
likely in the dense plasma in cores of giant elliptical galaxies be- 
cause, due to the high density, the equilibration time is short. If 
the higher-than-expected line ratios were due to incomplete atomic 
physics in the current spectral models, then they should not change 



as a function of radius. This can be directly verified in NGC 4636, 
which has the best combination of X-ray flux, temperature, and ex- 
posure. The data for this galaxy allow us to extract spectra and 
determine the {lui.os + hn.io)l hxiai line ratios in the combined 
spectrum from two regions surrounding the core at 15"-150" from 
the centre of the galaxy, in the crossdispersion direction of the RGS 
(see Fig. |2j. We follow here the same fitting procedure as for the 
cores of the galaxies. The best fitting results are shown in the third 
column of Table|2] The off-centre line ratios are consistent with the 
15.01 A line being optically thin in this region. The difference in 
the spectra of the core and surrounding regions of NGC 4636 can 
be seen clearly in Fig. [5] While the 15 A line is suppressed in the 
centre, it may be slightly enhanced outside the dense core. 

The best fit line broadening scale factor i- indicates that in 
NGC 4636 the distribution of the line producing ions is more 
peaked than the surface brightness, suggesting that the distribution 
of metals is centrally peaked. In the other systems this effect is 
smaller and less significant. 



5 MODELLING OF RESONANT SCATTERING IN 
NGC 4636 

To gain a better understanding of the observed line ratios, we sim- 
ulate the effects of resonant scattering on the radial intensity pro- 
files of the strongest X-ray emission lines in NGC 4636. For this 
we use the deprojected temperature and density profiles measured 
with Chandra. 

In our deprojection an alysis, we follow the approach described 
by IChurazov et all ( l2008h . We assume spherical symmetry, but 
make no specific assumption about the form of the underlying grav- 
itational potential. For a given surface brightness profile in n.^ an- 
nuli, we choose a set of {n^ ^ n^) spherical shells with the in- 
ner radii r(i), i = 1, . . . The gas emissivity £ is assumed to 
be uniform inside each shell, except for the outermost shell, where 
the gas emissivity is assumed to decline as a power law of radius: 
fi = fiout where Pom is a parameter. The deprojection process 
is a simple least square solution to determine the set of emissiv- 
ities in the set of shells (along with the emissivity normalization 
£out of the outer layers) that provides the best description of the ob- 
served surface brightness. The emissivity of each shell can then be 
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Table 3. Oscillator strengths and optical depths for the strongest X-ray lines in the spectrum of NGC 4636 for Mach numbers 0.0, 0.25, 0.5, 0.75. The expected 
suppression due to resonant scattering (///o)circ gives the value integrated within a radius of 2 kpc from the centre of the galaxy. The approximate values 
(///o)rgs were obtained by integrating within an "effective extraction region" which is 0.5' wide and 3' long. The metallicity is assumed to be constant with 
the radius. 
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Figure 6. Observed radial profiles of the electron density and temperature 
used to model resonant scattering in NGC 4636. The blue circles and red 
squares indicate the deprojected profiles determined from Chandra data 
with metallicity as a free parameter in the fit and with metallicity fixed 
to 0.68 Solar, respectively. The magenta triangles show the projected radial 
temperature profile with the metallicity as a free parameter in the fit. The 
solid lines show the parametrization of the gas density and gas temperature 
profiles given in eqs.|3]and|4] 



evaluated as an explicit linear combination of the observed quanti- 
ties. Since the whole procedure is linear, the errors in the observed 
quantities can be propagated straightforwardly. With our definition 
of rj (see subsection |3.2t . the projection matrix does not depend on 
energy. Therefore, the deprojection in any energy band is precisely 
that used to deproject the surface brightness. Thus, we can accu- 
mulate a set of spectra (corrected for background and readout) for 
each of the annuli, and apply the deprojection to determine the 
emissivities of each shell in each of the ACIS energy chan nels. 

The Chandra data were modeled using XS PEC V12 jAmaud 

Il996h and the APEC model dSmith et al.ll200lh . The gas temper- 
ature and normalization were free parameters in the model. The 
heavy element abundance was either a free parameter in the fit or 
fixed to 0.68 Solar (as determined from a spectral fit to Chandra 
data extracted within the radius of 1 ' from the centre of NGC 4636). 
The metallicity profile in NGC 4636 has large uncertainties, be- 



cause the metal abundance distribution in the deprojected Chandra 
spectra has a large scatter as a result of noise enhancement dur- 
ing deprojection coupled with the limited spectral resolution of the 
CCDs. 

The electron density profile was derived from the spectral nor- 
malization, fixing the proton to electron ratio to 0.83. Fig.|6]shows 
the deprojected electron density (lower panel) and the temperature 
(upper panel) as a function of distance from the centre of the galaxy. 
The circles and squares show the deprojected profiles for metal- 
licity as a free parameter and for metallicity fixed to 0.68 Solar, 
respectively. The triangles show the projected temperature profile 
with the metallicity as a free parameter. 

Based on the observed radial distributions of electron density 
We and temperature kT^ in NGC 4636, we adopt the following ap- 
proximate forms for the deprojected density and temperature pro- 
files: 



We = 1.2 X 10" 



and 



1 + 



[—1 
V0.25'/ 



(3) 



,r = 0.56i±ii^:^keV. (4) 

We consider two possible radial behaviors for metallicity: a 
simple constant abundance of 0.68 Solar; and a centrally peaked 
abundance distribution, which is more consistent with the best fit 
XMM-Newton RGS line profile {s = 0.40 ± 0.04). In the latter case 
we approximate the metallicity profile with the following function: 



2 + (r/0.8')' 

: 0.95 — ^ - 0.4. 

1 + (r/0.8')' 



(5) 



We calculate the optical depth of the 15 A line from the centre 
of the galaxy to infinity, t = J/ijCrodr, where is the ion concen- 
tration and CTo is the cross section at the line centre, which for a 
given ion is 



^JnhrgCf 



where the Doppler width is given by 

1/2 



+ 



tuib 
r-2 



(6) 



(7) 



In these equations re is the classical electron radius, / is the oscil- 
lator strength of a given atomic transition, Eq is the rest energy of a 
given line, A is the atomic mass of the corresponding element, njp is 
the proton mass, c is the speed of light, and V,urb is the characteristic 
velocity of isotropic turbulence. VtuA = V2ViD,iurb5 where is 
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the velocity dispersion in the line of sight due to turbulence. Using 
the adiabatic sound speed, c, = ^jkT /firrip, the expression for 
the broadening can be rewritten as 



A£d = Eo 



2kT, y 1 



1/2 



(8) 



where ^ = 0.6 is the mean particle mass, y is the adiabatic in- 
dex, which for ideal monatomic gas is 5/3, and M = Kmrb/t's is 
the corresponding Mach number. The line energies and oscillator 
strengths were taken from ATOMDeH and the NIST Atomic Spec- 
tra Databas^ 

The resonant scattering effect has been calculated usin g 
Monte-Carlo simulations, as described bv lChurazov et alj j2004) . 
We model the hot halo as a set of spherical shells and obtain 
line emissivities f or each shell using the APEC plasma model 
jSmith et alj|200ll) . We account for scattering by assuming a com- 
plete energy redistribution and dipole scattering phase matrix. The 
line of Fe xvn has a pure dipole scattering phase matrix. In Fig. |7] 
we show the optical depth of the 15.01 A Fexvii line as a func- 
tion of radius, for various turbulent velocities, assuming the peaked 
Fe abundance distribution. For higher turbulent velocities, the op- 
tical depth becomes smaller. In Fig. [8] we show the ratio of the 
15.01 A line intensity calculated including the effect of resonant 
scattering to the line intensity without accounting for resonant scat- 
tering for various turbulent velocities, both for constant and peaked 
iron abundance distributions. The figure clearly shows the effect of 
the resonant scattering: in the core of the galaxy the intensity of 
the line is suppressed, whereas in the surrounding regions the in- 
tensity of the line rises. The effect is more prominent in the case 
of the centrally peaked Fe abundance distribution, where the cen- 
tral suppression is larger and the enhancement in the outer part is 
stronger. 

We compare the models of resonant scattering with the ob- 
servations and determine the systematic uncertainties involved by 
simulating spectra corresponding to each of these models, which 
are then fit in the same way as the observed data. To this end, we 
multiply the surface brightness profile of NGC 4636 with the the- 
oretical suppression profiles for the 15.01 A line (Fig.[8]l to deter- 
mine the predicted RGS line profile for the four considered char- 
acteristic turbulent velocities, for both assumed abundance distri- 
butions. We simulate spectra with photon statistics comparable to 
that of the actual observation, using the best fit values for the core 
of NGC 4636 from Table [2] as input parameters, and convolving 
the 15.01 A line through these predicted line profiles. By fitting the 
data simulated with no turbulence, we obtain (///o)i5oia = 0-69 
and 0.71 for peaked and flat abundance profiles, respectively. The 
ratios of (///o)i5oia isotropic turbulent velocity of M = 0.25 
are 0.76 and 0.78, and for M = 0.50 they are 0.80 and 0.82, for 
peaked and flat abundance profiles, respectively. The systematic 
uncertainty on these values due to our lack of knowledge about the 
actual Fe abundance distribution is ±0.02. We add this uncertainty 
in quadrature to the uncertainties associated with our modeling of 
the line broadening, which is at most 0.03, resulting in a total un- 
certainty of 0.04. Thus we conclude that the turbulent velocities in 
NGC 4636 are relatively small and that isotropic turbulence with 
a characteristic velocity of M > 0.25 can be ruled out at the 90% 
confidence level; turbulent gas motions at M > 0.5 can be ruled 
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Radius, kpc 



Figure 7. Optical depth of the 15.01 A line calculated from the given radius 
to infinity, for isotropic turbulent velocities corresponding to Mach numbers 
0.0, 0.25, 0.5, and 0.75. A centrally peaked Fe abundance distribution is 
assumed. 



out at the 95% confidence level. (The sound speed in NGC 4636 is 
c, ~ 400kms-i.) 

In Table |3] we list the strongest lines observed in the spec- 
trum of NGC 4636, their oscillator strengths and optical depths for 
different Mach numbers. Clearly, Fexvii at 15.01 A is the most 
optically thick line in the spectrum, but Oviii ha and Fexviii at 
14.2 A also have relatively large optical depths. For M = 0, the 
depth of the Fexvii line at 17.05 A line is t = 0.5, meaning that 
our assumption that the 17.05 A and 17.10 A line blend is op- 
tically thin is not completely correct, and our absented values of 
{I/Io)is 01 A biased slightly high, which further increases the sig- 
nificance of our upper limits on the turbulent velocities. In Table[3] 
we also list the predicted suppression for the strongest emission 
lines, (///o)circ = /^d^ / /^dA, integrated within a circular region 
of radius 2 kpc, and within the RGS extraction region (///o)rgs- 
The suppression for the RGS is determined in a simplified way, 
integrating within our 0.5' wide, and for line emission effectively 
about 3' long, RGS extraction region. These predicted ratios in Ta- 
ble [3] were determined assuming a flat abundance distribution. By 
adopting an optically thin plasma in the analysis, the best fitting O 
abundance within the 2 kpc radius can be underestimated by ~10- 
25%, and the Fe abundance by -10-20%. 



6 DISCUSSION AND CONCLUSIONS 

High-resolution spectra obtained with XMM-Newton RGS reveal 
that the Fe xvii line at 15.01 A in the cores of the elliptical galaxies 
NGC 4636, NGC 1404, NGC 5813, and NGC 4472 is suppressed 
by resonant scattering. The effects of resonant scattering can be 
investigated in detail for NGC 4636. The comparison of the mea- 
sured suppression of the 15.01 A Fexvii line in the spectrum of 
NGC 4636 with the simulated effects of resonant scattering reveals 
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Radius, kpc 

Figure 8. Radial profiles of the ratio of the 15.01 A line intensity with and 
without the effect of resonant scattering, for isotropic turbulent velocities 
con'esponding to Mach numbers 0.0, 0.25, 0.5, and 0.75. The upper panel 
shows the profiles for a constant metallicity, and the lower panel for a cen- 
trally peaked Fe abundance distribution. 



that the characteristic velocity of isotropic turbulence in the core of 
the galaxy is smaller than 0.25 and 0.5 of the sound speed at the 
~90% and 95% confidence levels, respectively. The energy density 
in turbulence is: 



^tuib 



(9) 



where p is the density of the plasma. The thermal energy of the 
plasma is given as: 



_ 3 p 



(10) 



By combining these equations we obtain the ratio of turbulent to 
thermal energy: 



2 



(11) 



For the 90% and 95% confidence level upper limits of M < 
0.25 and 0.5 on the characteristic velocity of isotropic turbulence 
in NGC 4636, the upper limit on the fraction of energy in tur- 
bulent motions is 5% and 20%, respectively. This is in agree- 
ment with the upper limits on the non-thermal pressure com- 
ponent in elliptical g alaxies obtained on larger spatial scales by 
IChurazov et alj j2008l) and with the properties of AGN induced tur- 
bulence in three-dimensional simulations, that use adaptive mesh 
hyd rodynamic code with a subgrid model of turbulence and mix- 
ing jScannapieco & Briiggenll20()^ . 

A small value for the average turbulent velocity dispersion 
was also infer red in th e initia l analysis of the RGS spectra for 
NGC 4636 bv lXuetal.1 ( |2002|) . They estimated that the turbulent 
velocit ies in the galaxy d o not strongly exceed 1/10 of the sound 
speed. ICaon et all j2000l) studied the kinematics of line emission 
nebulae in NGC 4636 and found an irregular radial velocity curve 
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Figure 9. The deprojected temperature profile of NGC 4649 determined 
using Chandra data. It shows an unusual central peak. 



which indicates turbulent motions in the core. The measured ra- 
dial velocities of the emission line nebulae are mostly within about 
~ 100 km s"' . The kinematics of emission li ne nebulae in the sample 
of elliptical galaxies studied bv .Caon et afl boOO) indicates chaotic 
gas motions with velocities of about M = 0.2-0.4. In cooling cores, 
optical emission l ine nebula e often spatially coincide with soft X- 
ray emission (e.g. ISparks et al. 2004; Fabian et al. 2006), which in- 
dicates a coupling, suggesting that the two gas phases should share 
the same turbulent velocities. 

The dense cores of NGC 4636 and NGC 5813 show strong 
evidence of recent violent interaction between the hot plasma and 
the AGN. This interaction is expected t o indu ce turbulence through 
the mechanical activity of the outflows IChurazov et alj200lh . Our 
results indicate that the typical velocities of the AGN induced tur- 
bulent motions in the hot plasma on sub-kpc spatial scales are less 
than 100-200 km s"'. Gas motions on larger spatial scales likely 
have larger velocities. Observations of the Perseus cluster suggest 
both turbule nt and laminar gas m otions with velocities as high as 
700 km s ' jFabian et alj|2003bl lal). The observed lack of resonant 
scattering in the 6.7 keV Fe-K line within the inner ~100 kpc 
( IChurazovet"ai]|2004h of Perseus confirms the presence of such 
strong gas motions. AGN induced gas motions with large range of 
velocities along our line of sight, perhaps even on spatial scales 
smaller than a 1 kpc, could explain the observed lack of resonant 
scattering in the 15.01 A Fexvii line in M 87 (the analysis of 
a 0.5' wide extraction region centre d on the core of M 87 gives 
Uahm + /u7.io)//u5.oi = 1-3 ± 0.3, iWerneretalj|2006) . We note 
that the central AGN in Perseus and in M 87 are more active in 
X-rays than the AGN in the target galaxies of this study. 

The dense cores of galaxies showing smaller disturbance by 
AGN might be expected to exhibit less turbulence than the heav- 
ily disturbed systems. Therefore, we might expect the 15.01 A line 
to be more suppressed in the most relaxed galaxies in our sam- 
ple. While there is a weak indication of the suppression increasing 
from the least to the most relaxed systems: from NGC 5813 and 
NGC 4636, to NGC 4472, the large errorbars on the line ratios and 
the small number of galaxies in the sample do not allow us to draw 
firm conclusions. 

The best fit line ratio in the apparently relaxed galaxy 
NGC 4649 shows no indication for resonant scattering in its core. 
The lack of observable optical depth effects in this galaxy can be 
explained by the unusual temperature profile of its gaseous halo. 
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which exhibits a sharp temperature peak in its centre (see Fig. |9] 
Allen et al. in prep.). Within the central 0.1 kpc of the galaxy, the 
temperature increases by about a factor of 2 to ~ 1.5 keV, at which 
no Fe XVII is present in the gas, implying that all the observed Fe xvii 
emission is from lower density regions surrounding the hot core. 
The high temp erature might be du e to shock heating by the cen- 
tral AGN (e.g. lShurkin et aklllOOSl) . None of the other galaxies in 
our sa mple shows such an unusual temperature profile ( I Allen et aU 
I2OO6I Allen et al. in prep.). 

Unfortunately, the high temperatures of cluster cores do not 
allow us to use Fexvii lines as a diagnostic tool to measure the 
effects of resonant scattering and thus place firm constraints on tur- 
bulent velocities in the intra-cluster medium. No other lines have 
a comparable sensitivity to the turbulent velocities. However, the 
dense cores of elliptical galaxies, are in terms of the dominant 
physical processes like cooling and AGN feedback, sufficiently 
similar to cooling cores of clusters of galaxies to argue that the 
turbulent velocities on small spatial scales, and thus the turbulent 
pressure support in these systems, are also likely to be low. This 
is encouraging for cosmolog ical studies based on X-ray observa- 
tions of galaxy clusters (e.g.lAllen et alllioosl : iMantz et allbOOSl : 
iHenrv et al.ll2008l : IVikhlinin et al.ll2008h . Simulations predict that 
turbulent gas motions in the intra-cluster medium can provide 5%- 
25% of the total pressure support through the virial region of clus- 
ters, potentially bia sing low hydrostat ic mass estimates, even in re- 
laxed systems (e.g. iNagai et alj[2007l) . However, these predictions 
depend sensitively on the gas physics assumed. Our results argue 
that such turbulence may be suppressed in nature, at least within 
the inner regions of such systems. For the merging unrelaxed Coma 
cluster, a ~10% lower lim i t on tu rbulent pressure support was ob- 
tained bv lSchuecker et al.l ( l2004h from a Fourier analysis of pres- 
sure maps. In the range 40-90 kpc, the turbulent spectrum was 
found to be well described by a projected Kolmogorov/Oboukhov- 
type model. 

Future, deeper (of the order of ~500 ks) observations of large 
nearby elliptical galaxies with the XMM-Newton RGS should pro- 
vide even tighter constraints on turbulent velocities from resonant 
scattering. X-ray calorimeters on future missions such as Astro- 
H and International X-ray Observatory (1X0) will enable direct 
measurements of the velocity-broadening of emissi on lines in the 
hot g as in galaxies, groups, and clusters of galaxies ( Sunv aev et al.l 
l2003i : Inogamov & Sunyaev 2003; Rebusco et al. 2008), thereby 
allowing an independent check on our results as well as a criti- 
cal extension to higher mass systems. However, we note that even 
the excellent (2-5 eV) spectral resolution of these instruments will 
not allow measurements of turbulent motions with velocities sig- 
nificantly smaller than 100-200 km s"'. Even though, the broad 
point spread function of the Astro-H mirrors will prohibit investi- 
gations of the elfects of resonant scattering on spatial scales smaller 
than presented in this paper, it will for the first time allow us to di- 
rectly measure the gas motions and spatially map them out to larger 
radii. The next planned X-ray observatory able to produce a real 
breakthrough by mapping the gas motions in groups and clusters of 
galaxies at high spatial resolution, will be IXO. 

An important consequence of resonant scattering is the bias it 
can introduce in measurements of radial profiles of abundances of 
chemical elements. As we show in Table[3] the integrated emission 
of several strong resonance lines is significantly suppressed within 
the inner 2 kpc of NGC 4636. By ignoring the effects of resonant 
scattering, the Fe and O abundance will be underestimated by 10- 
20 % within this region. Thi s conclusion is in line with the results 
of ISanders &Fabian| ( |2006h . However, while resonant scattering 



makes abundances appear to be 10-20% lower than the real value, it 
cannot explain the much stronger abundance dips seen in th e cores 
of several clusters of galaxies (e.g. lSanders & Fabianll2002h . 
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